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ABSTRACT 

We utilize existing imaging and spectroscopic data for the galaxy clusters MS2137-23 and Abell 
383 to present improved measures of the distribution of dark and baryonic material in the clusters' 
central regions. Our method, based on the combination of gravitational lensing and dynamical data, 
is uniquely capable of separating the distribution of dark and baryonic components at scales below 100 
kpc. Our mass models include pseudo-elliptical generalized NEW profiles for constraining the inner 
dark matter slope, and our lens modeling takes into account both the ellipticity and substructure 
associated with cluster galaxies as necessary in order to account for the detailed properties of multiply- 
imaged sources revealed in Hubble Space Telescope images. We find a variety of strong lensing 
models fit the available data, including some with dark matter profiles as steep as expected from 
recent simulations. However, when combined with stellar velocity dispersion data for the brightest 
member, shallower inner slopes than predicted by numerical simulations are preferred, in general 
agreement with our earlier work in these clusters. For Abell 383, the preferred shallow inner slopes 
are statistically a good fit only when the multiple image position uncertainties associated with our 
lens model are assumed to be O'.'S, to account for unknown substructure. No statistically satisfactory 
fit was obtained matching both the multiple image lensing data and the velocity dispersion profile of 
the brightest cluster galaxy in MS2137-23. This suggests that the mass model we are using, which 
comprises a pseudo-elliptical generalized NEW profile and a brightest cluster galaxy component may 
inadequately represent the inner cluster regions. This may plausibly arise due to halo triaxiality or by 
the gravitational interaction of baryons and dark matter in cluster cores. The intriguing results for 
Abell 383 and MS2137-23 emphasize the need for a larger sample of clusters with radial arcs. However, 
the progress made via this detailed study highlights the key role that complementary observations of 
lensed features and stellar dynamics offer in understanding the interaction between dark and baryonic 
matter on non-linear scales in the central regions of clusters. 
Subject headings: gravitational lensing - galaxies:formation - dark matter 



1. INTRODUCTION 

Cold dark matter (CDM) simulations (both with and 
without the inclusion of baryonic physics) are a crucial 
tool and proving ground for understanding the physics 
of the universe on nonlinear scales. One of the most 
active aspects of research in this area concerns the form 
of the dark matter density profile. Key questions raised 
in recent years include: Is there a universal dark matter 
density profile that spans a wide range of halo masses? 
What is the form of this profile and how uniform is it 
from one halo to another? To what extent do baryons 
modify the dark matter distribution? 

Drawing o n a suite of N-body simulations, 
iNavarro et al.l ()1997l ) originally proposed that the 
dark matter density profiles in halos ranging in size 
from those hosting dwarf galaxies to those with galaxy 
clusters have a universal form. This 3-D density distri- 
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bution, termed the NEW profile, follows pdm oc 
within some scale radius, r^c, and falls off as poM oc r^^ 
beyond. Subsequent simulations indicated that the 
inner density profile could be yet steep er - poM r^^-^ 
(jMoore et all Il998t iGhigna et all \200(t ) . As computing 
power increases and numerical techniques improve, it is 
now unclear whether the inner dark matter distribution 
converges to a power law form rather than becom- 
ing progressi vely s hallower in slope at smaller radii 
(Power et aD 120031 : INavarro et all 120041 : iDiemand et~all 
12004 12005ir 

For comparisons with data, such simulations need to 
account for the presence of baryons. This is particu- 
larly the case in the cores of rich clusters. Although 
baryons represent only a small fraction of the overall clus- 
ter mass, they may be crucially important on scales com- 
parable to the extent of typical brightest cluster galax- 
ies. Much work is being done to underst and the likely 
interactions between baryon s and DM (iGnedin et al.l 
12004 iNagai fc Kravtsovll2005l : iFaltenbacher et al.ll2005| ). 
These simulations will provide refined predictions of the 
relative distributions of baryons and DM. 

This paper is a further step in a series which aims 
to present an observational analog to progress de- 
scribed above in the numerical simulations. At each 
stage it is desirable to confront numerical predictions 
with observations. Whereas some workers have made 
good progress in constraining the total density profile 
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(e.g. iBroadhurst et al.l ()2005bD ). in order to address the 
relevance of the numerical simulations we consider it im- 
portant to develop and test techniques capable of sep- 
arating the distribution s of d ark and baryonic compo - 
nents (e.g. Sand ct al. (20021: Zappacosta ct all (|2006f ): 
iBiviano fc Salu cci (2006): Mah davi et a l. ( 20071 )). 

This paper pr esents a ref i ned v ersion of the method 
fir st proposed bvlSand et al.l ()2002[ ). exploited more fully 
in lSand et all ([200^ (hereafter S04). S04 sought to com- 
bine constraints from the velocity dispersion profile of 
a central brightest cluster galaxy (BCG) with a strong 
gravitational lensing analysis in six carefully selected 
galaxy clusters in order to separate the baryonic and 
dark matter distributions. S04 carefully selected clusters 
to have simple, apparently 'relaxed' gravitational poten- 
tials in order to match broadly the 'equilibrium' status 
of the simulated dark m atter halos originally analyzed 
bv lNavarro et al.l (|1997f ) and subsequent simulators. For 
example, Abell 383 and MS2137-23 have almost circular 
BCGs (6/a=0.90 and 0.83 respectively), require a single 
cluster dark matter halo to fit the strong lensing con- 
straints (in contrast to the more typical clusters that 
require a multi-modal dark matter morphology - Smith 
et al. 2005), have previously published lens models with 
a relatively round dark matter halo {b/a— 0.88 and 0.78 
respectively - Smith et al. 2001; Gavazzi 2005), and dis- 
play no evidence for dynamical disturbance in the X-ray 
morphology of the clusters (Smith et al. 2005; Schmidt 
& Allen 2006). 

The merit of the approach resides in combining two 
techniques that collectively probe scales from the inner 
~10 kpc (using the BCG kinematics) to the ~100 kpc 
scales typical of strong lensing. Whereas three of the 
clusters contained tangential arcs, constraining the total 
enclosed mass within the Einstein radius, three contained 
both radial and tangential gravitational arcs, the former 
providing additional constraints on the derivative of the 
total mass profile. In their analysis, S04 found the gradi- 
ent of the inner dark matter density distribution varied 
considerably from cluster to cluster, with a mean value 
substantially flatter than that predicted in the numerical 
simulations. 

S04 adopted a number of assumptions in their anal- 
ysis whose effect on the derived mass profiles were dis- 
cussed at the time. The most important of these included 
ignoring cluster substructure and adopting spherically- 
symmetric mass distributions centered on the BCG. The 
simplifying assumptions were considered sources of sys- 
tematic uncertainties, of order 0.2 on the inner slope. 
Although the six clusters studied by S04 were carefully 
chosen to be smooth and round, several workers at- 
tributed the discrepancy between the final results and 
those of the simulations as likely to arise from these sim- 
plifying assumptions feartclm ann fc: M ene ghettil 120041 : 
iDaial fc Keetonll2003l : lMeneghetti et al.ll2005f ). 

The goal of this paper is to refine the data analysis for 
two of the clusters (MS2137-23 and AbeU 383) originally 
introduced by S04 using fully 2-D strong gravitational 
lensing models, thus avoiding any assumptions about 
substructure or spherical symmetry. The lensing mod- 
els are ba sed on an improved version of t he LENSTOOL 
program (Kneib 1993; Kneib et al. 1996; see Appendix; 
http://www.oamp.fr/cosmology/lenstool/). A major de- 
velopment is the implementation, in the code, of a 



pseudo-elliptical parameterization for the NEW mass 
profile, i.e. a generalization of those seen in CDM simu- 
lations, viz: 
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where the asymptotic DM inner slope is f3. This for- 
malism allows us to overcome an important limitation of 
previous work and takes into account the ellipticity of 
the DM halo and the presence of galaxy-scale subhalos. 
Furthermore the 2-D lensing model fully exploits the nu- 
merous multiply-imaged lensing constraints available for 
MS2137-23 and AbeU 383. 

The combination of gravitational lensing and stellar 
dynamics is the most powerful way to separate baryons 
and dark matter in the inner regions of clusters. How- 
ever, it is important to keep a few caveats in mind. 
Galaxy clusters are structurally heterogeneous objects 
that are possibly not well-represented by simple parame- 
terized mass models. To gain a full picture of their mass 
distribution and the relative contribution of their major 
mass components will ultimately require a variety of mea- 
surements applied simultaneously across a range of radii. 
Steps in this direction are already being made with the 
com bined use of strong and weak gravit ational lensing 
(e.g. lLimousin et all (|2006l ): lBradac et all ( ,2006,) ), which 
may be able to benefit furt her from informat i on pro vided 
from X-ray analyses fe .g. ISchmidt fc AllenI (l2006f) ) and 
kinematic studies (e.g. iLokas fc MamonI ( 20031 )'). A re- 
cent analysis has synthesized weak-lensing. X-ray and 
Sunyaev-Zeldovich observations in the cluster Abell 478 
- similar cross-disciplinary work will lend further insights 
into the mass distribution of clusters (jMahdavi et al.l 
|2003). 

Of equal importance are mass models with an appro- 
priate amount of flexibility and sophistication. For in- 
stance, incorporating models that take into account the 
interaction of baryons and dark matter may shed light 
into the halo formation process and provide more ac- 
curate representations of dark matter structure. Halo 
triaxiality, multiple structures along the line of sight and 
other geometric effects will also be important to charac- 
terize. At the moment, incorporating these complexities 
and securing good parameter estimates is computation- 
ally expensive even with sophisticated techniques such as 
the Markov-Chain Monte Carlo method. 

Numerical simulation results are often presented as the 
average profile found in the suite of calculations per- 
formed. Instead, the distribution of inner slopes would 
be a more useful quantity for comparison with individual 
cluster observations. Also, comparisons between simula- 
tions and observations would be simplified if projected 
density profiles of simulated halos along multiple lines of 
sight were to be made available. These issues should be 
resolvable as large samples of observed mass profiles are 
obtained. 

For the reasons above, comparing observational re- 
sults with numerical simulations is nontrivial. The ob- 
servational task should be regarded as one of develop- 
ing mass modeling techniques of increasing sophistica- 
tion that separate dark and baryonic matter, so as to 
provide the most stringent constraints to high resolution 
simulations which include baryons as they also increase 
in sophistication. The combination of stellar dynamics 
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and strong lensing is the first crucial step in this process. 
Its diagnostic power will be further enhanced by includ- 
ing other major mass components (i.e. the hot gas of 
the intracluster medium or the stellar contribution from 
galaxies) out to larger radii. 

A plan of the paper follows. In § [2] we explain the 
methodology used to model the cluster density profile by 
combining strong lensing with the BCG velocity disper- 
sion profile. In §[3] we focus on translating observational 
measurements into strong lensing input parameters. This 
section includes the final strong lensing interpretation of 
MS2137-23 and Abell 383. In §|4]we present the results of 
our combined lensing and dynamical analysis. In §[5] we 
discuss further systematic effects, limitations and degen- 
eracies that our technique is susceptible to - with an eye 
towards future refinements. Finally, in §[B]we summarize 
and discuss our conclusions. Throughout this paper, we 
adopt r as the radial coordinate in 3-D space and R as 
the radial coordinate in 2-D projected space. When nec- 
essary, we assume Ho=65 km s~^Mpc~^, 17m=0.3, and 
f)A=0.7. 

2. METHODS 

The intent of this work is to use the full 2D informa- 
tion provided by the deep Hubble Space Telescope (HST) 
imaging in two strong lensing clusters (MS2137-23 and 
Abell 383) in conjunction with the BCG stellar velocity 
dispersion profile in order to constrain the distribution 
of baryonic and dark matter. These two clusters were 
selected for further study from the larger sample pre- 
sented by S04 because, of the three systems with both 
radial and tangential gravitational arcs, these two pre- 
sented the shallowest DM inner slopes. 

2.1. Lens Modeling 

We use the updated LENSTOOL ray-tracing code to 
construct models of the cluster mass distribution. Our 
implementation of the mass profiles is identical to that of 
iGolse fc KneibI ()2002D . with the exception that we have 
generalized their pseudo-elliptical parameterization to in- 
clude ones with arbitrary inner logarith mic slopes. For 
the de tails, the reader is referred to both lGolse fc KneTbl 
()2002l ). and the Appendix. Here we briefly explain the 
lens modeling process and parameterization of the cluster 
mass model. 

Identifying mass model components and multiple- 
image candidates is an iterative process. Initially, multi- 
ple images are spectroscopically confirmed systems with 
counter images identified by visual inspection and with 
the aid of preliminary lens models, taking into account 
that gravitational lensing conserves surface brightness. 
Multiple images without spectroscopic confirmation were 
used in the case of Abell 383, since these additional con- 
straints helped clarify the role that galaxy perturber #1 
(Table [S]) played in the central regions of the cluster (see 
§ I3.3p . If the location of a counter image is tentative, 
especially if there are several possibilities or an interven- 
ing cluster galaxy confuses the situation, the system is 
not included in deriving the mass model. § 13.21 and § 13.31 
present a detailed description of the final multiple image 
list adopted. 

Once the multiple images are determined, the cluster 
mass model is refined and perturber galaxy properties 
are fixed. In general, a lens mass model will have both 



cluster and galaxy scale mass components. The cluster 
scale mass component represents the DM associated with 
the cluster as a whole plus the hot gas in the intracluster 
medium. In the limit that the cluster DM halo is spheri- 
cal (see Eqn IAGp , its density profile has the form of Eqn[TJ 
In the adopted parameterization, the DM halo also has 
a position angle (9) and associated pseudo-ellipticity (e) 
(see EqnES&ES])- 

Galaxy scale mass components are necessary to ac- 
count for perturbations to the cluster potential that seem 
plausible based on the HST imaging and are demanded 
by the observed multiple image positions. These compo- 
nents are described by pseu do- isothermal e lliptic al mass 
distributions (PIEMD; Kass iola fc Kovneil (|1993[ )V Each 
PIEMD mass component is parametrized by its position 
(xc, j/c), ellipticity (e), position angle (9), core radius 
(fcorc), cut-off radius (rcut) and central velocity disper- 
sion (cto)- The projected mass density, E, is given by: 

(2) 

where p'^—[(x— a;c)/(l-l-e)]^-|-[(?;— yc)/(l— e)]^ and the 
ellipticity of the lens is defined as e={a—b)/{a+b) ^. The 
total mass of the PIEMD is thus 3/27r(T§rcut/G. In order 
to relate Equation [2] to the observed surface brightness 
of the BCG in particular, we take E = (M*/L)/, where 
M^/L is the stellar mass to light ratio and / is the surface 
brightness, and find the following relation 

M,/L = 1.50^ aoVeut /{GL) (3) 

where L is the total luminosity of the BCG. The M^,/L 
of the central BCG will be used as a free parameter in 
our mass modeling analysis. Further details and prop- 
erties of the trunca t ed PI EM D model can be found in 
iNataraian fc KneibI (| 19971 ) and 'Limo usin et al.l (|2005D . 

The relevant parameters of the perturber galaxies (po- 
sition, ellipticity, core radius, cutoff radius and position 
angle) are assumed to be those provided via examina- 
tion of the HST imaging (see !^ l3.1l for details). Only the 
central stellar velocity dispersion, CToj is determined by 
optimization. At a particular stage in the process, the 
predicted multiple image positions are compared with 
those observed, and a value calculated (see § 12. 3p . 
The iteration stops when a minimum is reached. The 
sole criteria for adding a perturber galaxy was whether 
or not it was necessary for the lens model to match the 
multiple image positions. If adding an additional per- 
turber did not alter our interim minimum i we did not 
include it in our subsequent analysis. 

Two special cases were encountered during the above 
procedure. First, one of the galaxy perturbers in 
Abell 383 required a larger cutoff radius (rent) than im- 
plied from the light distribution. As is described in § 13. 3i 
this concentration is necessary to account for several of 
the multiple image positions. For this mass concentra- 
tion, not only do we determine the optimum tJo param- 
eter, but also the cutoff radius {rcut)- 

The other special case concerns the BCG in both 
galaxy clusters. These are assumed to be coincident 

* This quantity should not be confused with the quite different 
definition used for the pseudo-eUiptical generalized NFW profile, 
see the Appendix. 
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with the center of the cluster DM halo - one justified by 
the co-location of the BCG and centra l X-ray isophotes 
(iSmith et al.ll2005l : iGavazzi et al.l[2003l) . The BCG mass 
distribution, represented by a PIEMD model, comprises 
only the stellar mass. In this case the HST imaging is 
used to fix the BCG ellipticity and position angle, but 
since the measured stellar velocity dispersion is to be 
used as a constraint on the cluster mass profile, we leave 
the central velocity dispersion parameter (and hence the 
stellar Eqn ^ free in the tensing analysis. As 

the Jaffe density profile is used for the BCG dynami- 
cal analysis (S04), the PIEMD core and cutoff radius 
which best match the Jaffe surface density are adopted. 
§ 13.11 discusses further the results of the surface bright- 
ness matching between the Jaffe and PIEMD models in 
the clusters. 

2.2. Incorporating the Dynamical Constraints 

Apart from the use of the pseudo-elliptical gNFW pro- 
file for the dark matter component, the observational 
data and analysis methods adopted here are identical to 
those used by S04. In that work, the observed veloc- 
ity dispersion profile of the BCG was interpreted via the 
spherical Jeans equation (see Appendix of S04) to as- 
sist in the decomposition of the dark and baryonic mass 
components. This portion of the calculated by 

comparing the expected velocity dispersion profile of the 
BCG (which depends on the mass enclosed at a given ra- 
dius and the relative contribution of dark and luminous 
matter) given a mass model with the observed velocity 
dispersion profile, taking into account the effects of see- 
ing and the long slit shape used for the observations. 
Ellipticity in the BCG and its dark matter halo can be 
ignored as its effect on the veloci t y disp ersion profile will 
be small (e.g. iKronawitter et all (|2000D ). 

The reader is referred to S04 for the observational de- 
tails pertaining to the velocity dispersion profile, the sur- 
face brightness profile of the BCG and the subsequent 
dynamical modeling to constrain the cluster DM inner 
slope. 

2.3. Statistical Methods 

A x^-estiniator is used to constrain the acceptable 
range of parameters compatible with the observational 
data. First we use the strong lens model to calculate the 
likelihood of the lensing constraints and then we com- 
bine it with a dynamical model to include the kinematic 
information into the likelihood. 

The first step is the strong lensing likelihood. Once 
the lensing interpretation is finalized and the perturbing 
galaxy parameters are fixed, the remaining free param- 
eters are constrained by calculating a lensing over a 
hypercube which encompasses the full range of accept- 
able models, modulo a prior placed on the dark mat- 
ter scale radius (see § 12. 4p . In Bayesian terms this 
corresponds to adopting a uniform prior. The lens- 
ing value is calculated in the source plane identi- 
cally to that of (jSmith et al.ll2005l ). For each multiply- 
imaged system, the source location for each noted image 
ixmodei,i,ymodei,i) IS Calculated using the lens equation. 
Since there should be only one source for each multiple 
image set (with N images), the difference between the 
source positions should be minimized, hence: 



2 X ^ i'^model^i ■^mode/,i+l) ~^{yraodel,i yraodel,i+l) 

Xpos— 2^ 

1=1 s 

erg is calculated by scaling the positional error associated 
with a muliple image knot, <t/, by the amplification of the 
source A so that aj — Aag. The following analysis will 
assume two different positional errors for the multiply 
imaged knots, using uncertainties of aj —Q'.'2 and (X'S, 
referred to hereafter as the 'fine' and 'coarse' analyses, 
respectively. The case for each is justified below. 

The finer 0'.'2 error bar corresponds to the uncertainty 
in the multiple image knot positions as defined by the 
resolution and pixel size of the IIST/WFPC2 images. 
Excellent strong lensing fits are achieved with the finer 
0'.'2 error bar {x^/dof ~ 1), so that technically there is 
no need for increased uncertainties. The uncertainty is 
dominated by the spatial extension of the multiple im- 
age knots employed and the ability to identify surface 
brightness peaks. 

In contrast to our ability to match the image positions 
down to the resolution of the IIST/WFPC2 images, re- 
cent combined strong and weak lensing analyses of Abell 
1689 h a ve been unable to do s o (e.g . Broadhurst et ajj 
(I2005a[ ): iHalkola et all (j2006D : iLimousin et al. (2006)). 
Although Abell 1689 is a more complex cluster than those 
studied here, it does have the most identified multiple im- 
ages of any other cluster to date, and so can probe the 
overall mass profile on smaller scales and with many more 
constraints. As espoused in § 1, real galaxy clusters are 
complex systems that are likely not easily parameterized 
by simple mass models, and as the number and density 
of mass probes increases the more refined and complete 
the mass model necessary to match the data. In our case, 
where we have relatively few mass profile constraints (at 
least in comparison to Abell 1689), we adopt a coarser 
0'.'5 positional error which allows us to account for com- 
plexities in the actual mass distribution of our clusters 
that our small number of mass probes are insensitive to. 
This, plus the fact that we carefully chose our perturb- 
ing galaxies such that a lensing jdof ~ 1 was found 
should account for reasonable situations where we have 
missed an interesting perturber galaxy. By adopting too 
small a multiple image position uncertainty, the region in 
parameter space explored may be overly confined (such 
that, for example, the observed BCG velocity dispersion 
profile cannot be reproduced). 

The strong lensing analysis is performed with 5 free pa- 
rameters, analogous to those adopted in S04. These are 
the dark matter inner logarithmic slope (/3) , the pseudo- 
ellipticity of the potential (e) , the amplitude of the DM 
halo (^c), the dark matter scale radius (r^c), and the 
mass-to- light ratio of the BCG (M*/L). We choose to 
place a uniform prior on the dark matter scale radius, 
(r^c) based on past mass profile analyses of these clusters 
and results from CDM simulations in order to reduced 
computation time (see § 12. 4p . In practice, to evaluate 
the Xpos at each point in the hypercube, the pseudo- 
ellipticity of the cluster dark matter halo is optimized 
while simply looping over the remaining free parameters. 

Once the strong lensing x^ values are computed, atten- 
tion is turned to the dynamical data. In contrast to the 
strong lensing model, the dynamical model is spherically 
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symmetric and follows that presented by S04, with the 
value being: 



X. 



N 
i=l 



{<Ji,obs-cri .mode 
A? 



(5) 



where is the uncertainty in the observed velocity dis- 
persion measurements. 

The lensing and velocity dispersion values are 
summed, allowing for standard marginalization of nui- 
sance parameters and the calculation of confidence re- 
gions. 

2.4. Dark matter scale radius prior 

As mentioned in the previous sections, in order to limit 
computation time, a prior was placed on the dark matter 
scale radius r^c- This is justified both on previous mass 
profile analyses of these clusters and the results of CDM 
simulations. 

An array of CDM simulations has provided information 
not only on the inner dark matter density profile, but on 
the expected value of the scale radius, r^c, and its intrin- 
sic scatter at the galaxy clu s ter scale Ce.g. |Bullock et al.l 
£2004: Do lag et al.ll2004[ ). For ex- 
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ample. Bullo ck et al.l (12001) found that dark matter ha- 
los the size of small galaxy clusters have scale radii 
between 240 and 550 kpc (68% CL). ITasitsiomi et al.l 
((2004}, using higher resolution simulations with fewer 
dark matter halos found r^c of 450±300 kpc. Finally, 
iDolag et al.l ( f2004 ) studied the DM concentrations of 
galaxy clusters in a KCDM cosmology and found a typ- 
ical range of scale radii between r^c of 150 and 400 kpc. 
These results represent a selection of the extensive nu- 
merical work being done on the concentration of dark 
matter halos. 

Previous combined strong and weak lensing analy- 
ses of MS2137-23 have provided approximate va l ues fo r 
the scale rad i us (iGavazzi et al.l 120031 : iGavazzil I2005D . 
iGavazzi et al.l (|2003D found a best fitting scale radius of 
~130 kpc (and hints that the scale radius may be as 
low as ~70 kpc from their weak lensing data) for their 
analy sis of MS2137-23. A more recent analysis (|Gavazzil 
l2005l ) found a best fitting radius of ~170 kpc. Simi- 
larly for Abell 383, a recent X-ray analysis found a best- 
fitting dark matter scale radius of ~ 130 kpc (Zhang et 
al. 2007). 

Taking these observational studies into account, we 
chose a uniform scale radius prior between 100 and 200 
kpc for MS2137-23 and Abell 383, both for simpHcity 
and to bracket the extant observational results which of- 
ten have constraints at larger radii (and thus constrain 
the scale radius better) than the current work. It is worth 
noting that the extant observations of these two clusters 
indicate a scale radius which is on the low end of that 
predicted from CDM simulations. For a fixed virial mass, 
a smaller scale radius indicates a higher concentration, 
c = ryir/rsc- This could be due to the effect of baryonic 
cooling, which could increase halo concentration (as well 
as inner slope perhaps). It has been suggested that those 
halos with the highest concentration (again for a fixed 
mass) are those that are the oldest and with the least sub- 
structure, providing more indirect evidenc e that we have 
chosen 'relaxed' galaxy clusters to study (jZentner et al.l 



l2005f ). We briefly explore the consequences of changing 
our assumed scale radius prior range in § 15.2.21 

3. APPLICATION TO DATA 

We now turn to the observational data for MS2 137-23 
and Abell 383 and describe our methods for analyzing 
these in the context of lensing input parameters. 

3.1. BCG and Perturber Galaxy Parameters 

In order to fix the position angle and cUipticity of the 
perturber galaxies and BCG components, the IRAF task 
ELLIPSE is used to estimate the surface brightness profile 
at typically the effective radius. The measured parame- 
ters are fixed in the lensing analysis. The galaxy position, 
core radius (rcorc) and cutoff radius (rcut) are each cho- 
sen to match those fitting the photometry (Table [3]). For 
perturbing galaxies, this leaves only the PIEMD param- 
eter velocity dispersion (ctq) which must be adjusted to 
match the multiple imaging constraints, as explained in 

§0 

For the BCG, following S04, it is preferable to use the 
Jaffe stellar density profile for the dynamical analysis 
since this function provides an analytic solution to the 
spherical Jeans equation. However, the PIEMD model 
implemented in LENSTOOL offers numerous advantages 
for the lensing analysis. To use the most advantageous 
model in each application, a correspondence is estab- 
lished between the two by fitting with a PIEMD model 
the Jaffe surface brightness fit presented by S04. An ap- 
propriate combination of the PIEMD rcoro and Tcut model 
parameters matches the Jaffe profile found by S04 with 
no significant residuals (PSF smearing was also taken 
into account). Table [T] lists the PIEMD parameters used 
for our lensing analysis, as well as the Jaffe profile pa- 
rameters used by S04. 

3.2. MS2137-23 Lens Model 

The strong lensing properties of MS 2137-23 have 
been st udied extensively by ma ny workers (iMellier et al.l 
biggs'; 'Miralda-Escude 1995^; jHammer et al.l Il997t 
Gavazzi et al. 2003; Gavazz i 2003). The most detailed 
model (jCavazzi et al.l l2003f ) used 26 multiply-imaged 
knots from two different background sources. The model 
adopted here is more conservative and based only on 
those multiple images confirmed via spectroscopy or 
suggested on the grounds of surface brightness or interim 
lens models. Despite having some multiple images in 



common with IGavazzi et air(|2003l ). we have retained 



our own no menclature. 

Following ISand et "aTl (|2002f ). the tangential and ra- 
dial arcs arise from separate sources, at z — 1.501 and 
z = 1.502, respectively. The multiple image interpre- 
tation is detailed in Figure [1] and Table [H There are 
two separate features (1 and 3) on the source giving rise 
to the tangential arc which is multiply- imaged four and 
three times respectively. It has not been possible to confi- 
dently locate the fourth image of feature 3, since it is ad- 
jacent to the perturbing galaxy. Also, it is expected that 
a fifth, central image would be associated with the giant 
tangential arc. Although the position of this fifth cen- 
tral i mage has been tentatively reported (IGavazzi et al.l 
I2OO3I ). we do not include it in our model because we 
were unable to clearly identify it due to BCG subtrac- 
tion residuals. Two images of the source giving rise to 
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TABLE 1 

Fixed Parameters in Abell 383 and MS2137-23 Lens Model 



Cluster 


Parameters 


(arcsec) 


Vc 

(arcsec) 


D/a 


71 

(deg) 


(kpc) 


O"0 
(kms ) 


^ 71? 

rent/ -tie 

(kpc) 


iVlozio ( 


Cluster-scale DM halo 


n n 
U.U 


n n 
U.U 


Free 


o.U 










BCGpiEMD 


0.0 


0.0 


0.83 


17.75 


5x10-'^ 


Free 


22.23 




BCGjaffc 


0.0 


0.0 










24.80 




Galaxy Perturber 


16.2 


-5.46 


0.66 


159.9 


0.05 


173.0 


4.8 


Abell 383 


Cluster-scale DM halo 


0.0 


0.0 


Free 


104.3 










BCGpiEMD 


0.0 


0.0 


0.90 


107.2 


3x10"'' 


Free 


25.96 




BCGjaffc 


0.0 


0.0 










46.75 




Perturber 1 


14.92 


-16.78 


0.804 


-20.9 


0.67 


230.0 


26.98 




Perturber 2 


9.15 


-1.92 


0.708 


10.3 


0.51 


140.0 


10.79 




Perturber 3 


0.17 


-24.26 


0.67 


65.2 


0.24 


124.8 


9.10 




Perturber 4 


-4.10 


-13.46 


0.645 


27.7 


0.17 


125.7 


2.19 



Note. — The position angle, 9 is measured from North towards East. The DM halo is par ameterized with the pseudo-gNFW profile. 
All other mass components are parameterized by a PIEMD model. Note that Re = 0.7(irjaffe l lJaffelll983l) . 



the radial arc were also identified. The mirror image of 
feature 2a nearer the center of the BCG could not be 
recovered, most likely because of residuals arising from 
the sub traction of the BCG . 

As in lGavazzi et al.l ()2003f ). only one perturbing galaxy 
is included in the lens model (see Table [1]). For this 
system, the best- fitting value occurred near ao—l73 
km s~^. 

In the initial modeling of MS2137-23, some experimen- 
tation was undertaken with different cluster DM elliptic- 
ities and position angles. While some variation in ellip- 
ticity is permitted by the lensing interpretation, a robust 
position angle offset was detected between the BCG and 
that of the DM ha l o of A O ~13 degrees, in agreement 
with iGavazzi et aTl (|2003f ). In the following, results are 
presented with the DM position angle fixed at 6=5 de- 
grees. This optimal position angle was determined dur- 
ing the initial lens modeling process by fixing all cluster 
mass parameters to values corresponding to a model with 
X^/dof ^ 1 and letting the DM position angle vary until 
a minimum was reached. As a consistency check, we 
repeated our calculations with a fixed DM position angle 
of 4.0 and 6.0 degrees. Varying the DM position angle 
had very little effect on our other parameter constraints, 
but results in a slightly larger overall (lensing -I- veloc- 
ity dispersion profile; A^^ <1) value. For this reason, we 
only present our DM position angle of 5 degree results. 

3.3. Abell 383 Lens Model 

Detailed lens mode l s for A bell 383 have been published 
in lSmith et al.l (|200lL 120051 ) . which we will largely adopt 
in this work. Multiple image sets 1 and 2 are based o n 
the in-depth lensing interpretation of lSmith et al.l (|2005D . 

The reader is referred to that work for a detailed de- 
scription of this radial and tangential gravitational arc. 
Multiple image sets 3, 4, 5 and 6 (for which there are 
no spectroscopic redshifts, but for which their distinc- 
tive morphology is reassuring) are included largely to 
constrain the properties of perturbing galaxies 1,3, and 
4 (see Fig [2] and Table [2]). Since these images have no 
spectroscopic confirmation, a redshift z ~3 was assumed; 
the mass model is very insensitive to the exact choice. 

The Abell 383 cluster mass model is more complex 
than that for MS2137-23, but only in the sense that 
there are more perturbing galaxies that must be put 
into the mass model to match the image positions. The 
bright cluster elliptical southwest of the BCG (#2 in 



ISmith et all (|2005[ ): Perturber #1 in this work) requires 
a DM halo more extended than the light, as mentioned 
in § 12.11 After some iteration, it was found that the pa- 
rameters of this important perturber could be fixed to 
those listed in Table [1] Multiple image sets 3, 4 and 5 
play a crucial role (see Table [2|) in constraining the per- 
turber. Although other perturbing galaxies were added, 
none has a comparable effect on the lensing x^ . Table [T] 
provides the full model parameter list. 

A slight (^--^3 degree) offset between the position angle 
of the BCG and the cluster DM halo was noted. We 
found the best-fitting DM position angle in the same way 
as in MS2137-23 (§|3l2]). The position angle of the DM 
halo was kept fixed but the ellipticity was left as a free 
parameter. As in MS2137-23, we also reran our analysis 
with a DM position angle of 103.3 and 105.3, but only 
present the results with a DM PA of 104.3. 

4. RESULTS 

We now analyze the refined 2D lens models of MS2137- 
23 and Abell 383 together with the velocity dispersion 
profiles presented in S04. We present our analysis with 
both a multiple image position uncertainty of 0'.'2 and 
0'.'5, described as the fine and coarse fits in § 12.31 The 
results are summarized in Table |3] and Figures [3] and IH 

4.1. MS2137-23 

Figure [3] and the discussion below summarizes the re- 
sults for the fine and coarse positional cases. One thing 
to note is the number of degrees of freedom involved, 
i.e. the difference between the number of constraints 
and the number of free parameters, in order to quantify 
the goodness of fit. The mass model has 5 free param- 
eters, as detailed in § 13.21 the DM inner slope /3, the 
DM pseudo-ellipticity e, the DM amplitude Sc, the BCG 
stellar M^/ L, and the dark matter scale radius Tsc which 
is allowed to vary in the 100-200 kpc range. Consider- 
ing Table [21 the multiple images provide 12 constraints, 
while the velocity dispersion data provides 8, giving a 
total of 20 data constraints. The resulting number of 
degrees of freedom is thus 15. 

4.1.1. The fine positional accuracy lensing case 

The two top panels of Figure [3] and the appropriate 
line in Table [3] encapsulate the results of the fine posi- 
tional analysis. DM inner slopes between (3 = 0.65 — 1.0 
lie within the 68% confidence limit (after marginalizing 




Fig. 1. — Multiple image interpretation of the cluster MS2137-23. The exact positions used arc shown in Table[2] Three sets of multiple 
images are identified, one with the radial arc system (2a & 2b) and two with the tangential arc system (labcd & 3abc). The perturbing 
galaxy is the elliptical SO next to the lensed feature lb. 



TABLE 2 

Multiple Image Interpretation of MS2137 and Abell 383 



Cluster 


Multiple Image 






Redshift 




ID 


(arcsec) 


(arcsec) 




MS2137 


la 


6.92 


-13.40 


1.501 




lb 


12.40 


-7.94 


1.501 




Ic 


0.07 


19.31 


1.501 




Id 


-11.57 


-7.49 


1.501 




2a 


3.96 


-5.51 


1.502 




2b 


-8.01 


22.10 


1.502 




3a 


5.16 


-14.68 


1.501 




3b 


0.11 


18.91 


1.501 




3c 


-12.30 


-6.74 


1.501 


Abell 383 


lA 


-1.74 


2.56 


1.0 




IB 


-1.03 


1.20 


1.0 




IC 


16.37 


-4.03 


1.0 




2A 


7.00 


-14.01 


1.0 




2B 


8.23 


-13.20 


1.0 




2C 


14.11 


-8.19 


1.0 




3A 


5.88 


-22.02 


3.0 




3B 


14.69 


-14.68 


3.0 




3C 


16.49 


-14.39 


3.0 




4A 


8.35 


-23.96 


3.0 




4B 


17.45 


-17.28 


3.0 




4C 


17.92 


-15.43 


3.0 




5A 


6.64 


-21.75 


3.0 




5B 


16.98 


-14.09 


3.0 




6A 


7.05 


-21.75 


3.0 




6B 


17.27 


-14.17 


3.0 



Note. — Multiple Image Interpretation. All image positions are with respect to the BCG center. 
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Multiple image interpretation of the cluster AbcU 383. The exact positions used are shown in Table [2] In all figures except 
we have subtracted cluster galaxies in order to more clearly see the multiple image features. 
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5 4 6 1 10 iOO 

M/L, K (kpo) 

Fig. 3. — The combined lensing+dynamics results for the cluster MS2137-23. The top row summarizes the results for the 0'.'2 lensing posi- 
tion uncertainty scenario while the bottom row encapsulates the 0'.'5 scenario. Top Left— Lensing+dynamics likelihood contours (68%, 95%, 
and 99%) in the M/L — /3 plane after marginalizing over the other free parameters with the 0'.'2 lensing multiple image uncertainty. Top 
Right- Best fitting velocity dispersion profile from the combined lensing-|-dynamics analysis with the 0''2 lensing multiple image uncer- 
tainty. No models could be found that fit both the lensing and observed velocity dispersion constraints. Bottom Left — Lensing+dynamics 
likelihood contours (68%, 95%, and 99%) with the 0"5 lensing multiple image uncertainty in the M/L — f3 plane after marginalizing or 
optimizing over the other free parameters. Bottom Right - Best fitting velocity dispersion profile from the combined lensing+dynamics 
analysis with the 0''5 lensing multiple image uncertainty. While the best-fitting model velocity dispersion is a better fit to the data than 
in the 0'.'2 lensing scenario, it still cannot reproduce the observed decline in the velocity dispersion profile in our highest radial bins — 
suggesting a problem with our current mass model parameterization. 



over all other free parameters), although the best fitting 
DM scale radius sits at the edge of the allowed prior 
(?'sc. best =200 kpc). A scale radius of 200 kpc is higher 
than that seen in previous lensing analyses of MS2137- 
23 that have used a sim i lar m a ss param e teriza tion to our 
own (e.g lGavazzi et all ()2003[ ): iGavazzil (|2005f) ). so there 
is no case to alter the prior. 

The parameter constraints are not particularly tight 
because the total ^54, larger than expected given 
the number of degrees of freedom. Such a value may 
indicate that the form of the mass profile used in the fit 
is inappropriate. Indeed, the model velocity dispersion 
profile is a poor match to that observed (Figure [3]). In 
fact, if the BCG velocity dispersion results are ignored, 
acceptable lens models {x^/dof ^1) can be recovered 
with a variety of inner DM slopes, scale radii and BCG 
stellar M/L, although these parameters have correlated 
values. We postpone discussion of the possible reasons 
for this mismatch until later. 

4.1.2. The coarse positional accuracy lensing case 

The two bottom panels of Figure [31 along with Table [3] 
summarize our results for the coarse positional analysis. 
DM inner slopes between f3 =0.4-0.75 lie within our 68% 
CL, and we again find DM scale radii at the high end of 
our prior range, as expected. The shift in the BCG M/L 
vs. DM inner slope contour fine positional case indicates 
that the increased parameter space in the lensing models 
has led to a slightly improved velocity dispersion profile 
(bottom right panel of Figure [3]). The overall —31 



is improved, although the probability for 15 degrees of 
freedom is less than 1%, assuming measurements are gov- 
erned by Gaussian statistics. Thus the model remains a 
marginal fit to the data. 

4.1.3. Comparison with Gavazzi 2005 

We n o w br iefly compare our results with those of 
IGavazzil ()2005l ). Gavazzi's analysis used a similar strong 
lensing model to our own, including what we consider 
to be somewhat less secure multiple images. However, 
he extended the analysis to larger scales including weak 
lensing data and incorporated the BCG velocity disper- 
sion profile presented by S04. Gavazzi adopted a strict 
NFW profile for the cluster DM halo and a Hernquist 
profile for the stellar component of the BCG. 

Despite these differences, Gavazzi's conclusions are 
very similar to those of the present paper. Models with 
NFW-like DM haloes (regardless of whether the inner 
slopes were varied) were uniformly poor fits to the obser- 
vational data. In particular, the falling velocity disper- 
sion profile observed at i? > 5kpc cannot be reproduced, 
despite experimenting with the effect of anisotropic or- 
bits in the stellar distribution. A major conclusion of 
Gavazzi's study is that halo triaxiality, an effect not typ- 
ically included, may play an important role in the central 
regions of galaxy clusters. We will return to this topic in 

4.2. Ahell 383 
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TABLE 3 
Acceptable Parameter Range 



Prior Setup Cluster Inner DM Slope e Af* /Ly 
^ (kpc) 



100 - 


200 kpc Tsc Prior 
















=0'.'2 


MS2137 


95+0-05 






200-*2 








Abell 383 


"S^+o-^ 
^■^''-o.os 


0.08«;«^ 


1 40000+^500 


100+28 


2 4+0-42 
^-^-0.42 




=0'.'5 


MS2137 


&+°'^^ 

U.O_g 2 


0.06«;«1 


44600+^^500 


200-31 


2 45+0-75 






Abell 383 


45+0-2 


0.06«:«? 


156000tf,500 


100+21 


2 34+1.02 



Note. — Best- fitting parameters and/or confidence limits for the difi'erent prior scenarios present in this paper. 



Our results for Abell 383 are shown in Figure|l]and Ta- 
ble [3] for both the fine and coarse positional uncertainty 
cases. We discuss each separately belo-w. Calculating the 
number of degrees of freedom in a similar -way to that 
done for MS2 137-23, we again have 5 free parameters in 
our mass model. Considering Table [21 multiple images 
provide 16 constraints, taking into account that those re- 
lated to multiple image sets 3,4,5 & 6 do not have known 
redshift information. The velocity dispersion data pro- 
vide 3 additional constraints. Thus, the resulting number 
of degrees of freedom is 14. 

4.2.1. The fine positional accuracy lensing case 

The top two panels of Figure 3] and the appropriate line 
in Table [3] summarize the results in this case. DM inner 
slopes between (3 = 0.5—0.7 lie within the 68% confidence 
limit of our analysis (after marginalizing over all other 
free parameters). The best fitting scale radius sits again 
at the edge of the allowed r^c prior range (rsc=100 kpc). 
An X-ray analysis of Abell 383, which was able to probe 
to higher radius than the current analysis, indicates that 
the DM scale radius is well above 100 kpc (Zhang et al. 
2007). For these reasons, and those discussed earlier, 
there is no case for adjusting the DM scale radius prior. 

The total x^=40.4, high given the 14 degrees of free- 
dom in the analysis. 

4.2.2. The coarse positional accuracy lensing case 

The two bottom panels of Figure [H along with Table [3] 
summarize the results for the coarse positional accuracy 
case. DM inner slopes between (3 =0.2 - 0.65 lie within 
our 68% CL, along with a best fitting DM scale radius 
of 100 kpc. Our parameter constraints encompass the 
values found in the fine accuracy case with no shift in 
parameter space (unlike the case for MS2137-23). This 
suggests that although we should expect a lower due 
to the increased uncertainties allowed, no significant im- 
provement to the best-fitting velocity dispersion profile 
should be expected. As we can see in the bottom right 
panel of Figure |4l the best fitting velocity dispersion pro- 
file is very similar to that obtained in the fine case. The 
total x^=22, acceptable given 14 degrees of freedom. 

5. DISCUSSION 

In the previous section we have presented the results of 
our analysis, which showed that a mass model comprising 
a stellar component for the BCG following a Jaffe pro- 
file together with a generalized NFW DM cluster halo is 
able to adequately reproduce the observations for Abell 
383 (albeit only for the coarse lensing positional accuracy 
scenario) but is unable to simultaneously reproduce the 



observed multiple image configuration and BCG veloc- 
ity dispersion profile for MS2137-23. In the case of Abell 
383, the inner DM profile is flatter than /3=1, supporting 
the earlier work of S04. This indicates that at least some 
galaxy clusters have inner DM slopes which are shallower 
than those seen in numerical simulations - but only if the 
mass parameterization used in the current work is reflec- 
tive of reality. Further work in this interesting cluster 
using other observational probes will further reflne the 
mass model, and determine if the generalized NFW DM 
form is a good fit to the cluster profile. 

In this section we discuss systematic uncertainties in 
our method and possible refinements that could be made 
to reconcile the mass model with the observations for 
MS2137-23. We hope that many of these suggestions 
will become important as cluster mass models improve 
and thus will present fruitful avenues of research. 

5.1. Systematic Errors 

We focus first on systematic errors associated partic- 
ularly with the troublesome stellar velocity dispersion 
profile for MS2137-23. Errors could conceivably arise 
from (i) significant non-Gaussianity in the absorption 
lines (which are fit by Gaussians) , (ii) uncertain measure- 
ment of the instrumental resolution used to calibrate the 
velocity dispersion scale, and (iii) template mismatch. 

Non-gaussianity introduces an error that we con- 
sider too small to significantly alter the goodness of fit 
(jCavazzi 2005) . The instrumental resolution of ESI (the 
Keck II inst rument used to mea sure the velocity disper- 
sion profile; ISheinis et al.l (|2002D 1 is ~30 km s~^; this is 
much smaller than the measured dispersion. Even if the 
instrumental resolution was in error by a factor of two, 
the systematic shift in a would only be 3 km s^^ (using 
Eq 3 in Treu et al. 1999). This would affect all measure- 
ments and not reverse the trend with radius. 

Concerning template mismatch, S04 estimated a pos- 
sible systematic shift of up to 15-20 km s~^. This could 
play a role especially as the signal to noise diminishes at 
large radii, where the discrepancy with the model profile 
is greatest. To test this hypothesis, we added 20 km s~-'^in 
quadrature to only those velocity dispersion data points 
in MS2137-23 at i? > 4 kpc and recalculated the best- 
fitting values, is reduced from 31 to 28.8, a modest 
reduction which fails to explain the poor fit. 

Although selectively increasing the error bars on those 
data points most discrepant with the model is somewhat 
contrived, our result does highlight the need for high S/N 
velocity dispersion measures out to large radii. A high 
quality velocity dispersion profile h as been measu red lo- 
cally for Abell 2199 to -20 kpc (|Kelson et all [2002.) . 
Interestingly, these high S/N measures display similar 
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5 4 6 1 10 iOO 

M/L, K (kpo) 

Fig. 4. — The combined lensing+dynamics results for the cluster Abell 383. The top row summarizes the results for the 0''2 lensing position 
uncertainty scenario while the bottom row encapsulates the O'.'S scenario. Top Left-Lensing+dynamics likelihood contours (68%, 95%, 
and 99%) in the Af/L — /3 plane with the 0'.'2 lensing multiple image uncertainty after marginalizing over the other free parameters. 
Top Right- Best fitting velocity dispersion profile from the combined lensing+dynamics analysis with the 0'.'2 lensing multiple image 
uncertainty. Bottom Left - Lensing+dynamics contours (68%, 95%, and 99%) in the Af/L — /3 plane with the 0'.'5 lensing multiple image 
uncertainty after marginalization over the other free parameters. Bottom Right - Best fitting velocity dispersion profile from the combined 
lensing+dynamics analysis with the 0''5 lensing multiple image uncertainty. The 0'.'5 lensing multiple image scenario provides a better 
overall fit to the observations, although we are limited by the relatively poor quality of the observed Abell 383 velocity dispersion profile. 



trends to those for MS2137-23 in the overlap regime, i.e. 
a slightly decreasing profile at i? < lOkpc. The dip wit- 
nessed in MS2137-23 is thus not a unique feature, al- 
though with deeper measurements we might expect to 
see a rise at larger radii as a result of the shallow DM 
profile. 

A final potential limitation in the dynamical analy- 
sis is t he assumption of orbital isotropy. Both S04 and 
iGavazz i (2005) explored the consequences of mild orbital 
anisotropy, concluding a possible offset of A/3 ~ 0.15 
might re sult. Even inclu ding orbital anisotropy into his 
analysis, iGavazzil (|2005[ ) was unable to fit the observed 
velocity dispersion profile. 

Since we determine our lensing values in the source 
plane, we checked to make sure that no extra images were 
seen after remapping our best-fit lensing -t- velocity dis- 
persion models back to the image plane. No unexpected 
images were found, although several images that were 
explicitly not used as constraints were found, such as the 
mirror image of radial arc image 2a in MS2137 and the 
complex of multiple images associated with 3abc, Sab, 
and 6ab in Abell 383 (see Figures [1] and [2|) . As discussed 
in § 13.21 and 13.31 some of these multiple images were 
not used as constraints because we could not confidently 
identify their position either due to galaxy subtraction 
residuals or blending with other possible multiple image 
systems. 

We finally comment on the uncertainties assigned to 
the multiple image systems for our lens models. We 
have presented two sets of results in this work; with as- 
signed image positional accuracies of o'/=0'.'2 and O'.'S. 
We find a variety of lens models are compatible with the 



iT/=0'.'2 case and only when the velocity dispersion data 
is included into the analysis does the data fail to be re- 
produced by the model. Certainly if we were to further 
increase the positional errors, at some point a good ve- 
locity dispersion fit could conceivably be obtained, but 
we will refrain from doing so in the present work. 

Increasing the positional uncertainties is only justified 
if there is evidence that there are significant missing com- 
ponents in the mass models. Further observations that 
can probe the mass distribution on fine scales to larger 
radii and higher quality models which can account for 
phenomena such as adiabatic contraction in the inner re- 
gions of galaxy clusters and triaxiality represent the best 
way to obtain a more precise picture of the cluster mass 
distribution. 

5.2. Improving the Mass Model 

We now turn our attention to possible inadequacies 
in the mass model. It is important to stress that the 
two diagnostics (lensing and dynamics) adopted in this 
study probe different scales. The lensing data tightly 
constrains the mass profile at and outside the radial arc 
(^--^20 kpc), while the velocity dispersion constrains the 
mass profile inside i? < 10 kpc. Since multiple images are 
numerous and their positions can be more precisely mea- 
sured than velocity dispersion ^, they carry more weight 
in the statistic than the kinematic points, producing 
a best overall fitting model (which is lensing dominated) 

^ The error on the astrometry with respect to the relevant scale, 
the Einstein Radius Se is much smaller than the relative error on 
velocity dispersion, i.e. 59/9^ << Scr/a 
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r„=400kpc, fixed 


MS2137-2353 
Moore 


NFW 









MS2137-a353 
rj^ = 50kpc, fixed Moore 



Abell 383 
rj^=50kpc, fixed Moore 



NFW 



M/L, M/L, 

Fig. 5. — Confidence contours (68%, 95%, and 99%) when we allow the dark matter scale radius to be fixed at values a factor of two 
beyond our observationally motivated prior. Top Row - Contours when we fix the dark matter scale radius to rsc=50 and 400 kpc in 
MS2137. Although the rsc=400 kpc scenario provides a relatively good fit to the data (x^ ~26), this value for the scale radius is much 
larger than that observed from weak lensing data. The rsc=50 kpc scenario is a significantly worse fit to the data, with ~39. Note 
that the DM inner slope is /3 < 1 in both scenarios. Bottom Row - Contours when we fix the dark matter scale radius to rsc=50 and 400 
kpc in A383. The large discrepancy in inner slope values obtained emphasize the need for a mass probe at larger radii. The best-fitting 
model for either fixed scale radius is significantly worse than the best-fitting rsc=WO kpc result (x^ ~26.5 and 31.3 for rsc=50 and 400 
kpc respectively). 



Abell 383 
r^„=400kpc, fixed Moore 




that is a relatively poor fit to the kinematic data. To im- 
prove the model, one must admit that either one of the 
two components of the modeling is incorrect, or that the 
functional form of the mass profile chosen to extrapolate 
the lensing information at the scales relevant for dynam- 
ics is insufficient. In this section we discuss several areas 
where the mass model presented in this paper could be 
improved. 

5.2.1. The Contribution of the Brightest Cluster Calaxy 

We might query the assumption of a Jaffe density pro- 
file for the BCG. This seems an unlikely avenue for im- 
provement given the Jaffe profile fits the observed BCG 
surface brightness profile re r narkab ly well (see Figure 2 
of S04) . Moreover, iGavazzil (j2005l ) utilized a Hernquist 
mass profile in his analysis o f MS213 7 -23, w hich also 
matches the observations, and IGavazzil ()2005D was like- 
wise unable to reproduce the observed S04 velocity dis- 
persion profile. 

We have additionally checked our assumptions by al- 
tering the PIEMD fit to the BCG surface brightness data 
so that it is matched not to the derived Jaffe profile fit 
to the BCG but directly to the HST surface brightness 
profile. With this setup, we found a Vcut value of 23.70 
kpc for MS2137 and 28.65 kpc for AbeU 383 (compare 
this with the numbers in Table [Ij . Redoing our analy- 
sis for the best-fitting r^c scenario only, our constraints 
on (3 for both AbeU 383 and MS2137 did not change by 
more than 0.05, and so it is not likely that our method 
for constraining the BCG mass contribution is the root 
cause of our inability to fit the data to a BCG -I- gNFW 
cluster DM halo mass model. 



Conceivably the BCG may not be coincident with the 
center of the cluster DM halo, as has been assumed 
throughout this work. It is often the case that small sub- 
arcsecond off sets between BCGs and clust er DM halos 
are ne cessary to fit lensing constraints (e.g. ISmith et al.l 
(120051) '). There is strong evidence that the BCG is nearly 
coincident with the general cluster DM halo in projec- 
tion fro m the stron g lensing work presented here and by 
others (iGavazzi et a l. 2003; Gavazzi 2005). However, an 
offset could be responsible for the flat to falling observed 
velocity dispersion profile if the BCG were actually in a 
less dark matter dominated portion of the cluster. An- 
other possibility is that there are multiple massive struc- 
tures along the line of sight, which would be probed by 
the strong lensing analysis, but not with the velocity dis- 
persion profile of the BCG. A comprehensive redshift sur- 
vey of MS2137-23 could provide further information on 
structures along the line of sight. 

5.2.2. The Advantage of a Mass Probe at Larger Radii 

With our presented data set, we have seen that it is 
difficult to constrain the DM scale radius, r^c because 
both of our mass probes are only effective within the 
central ^^100 kpc of the clusters - within the typical DM 
scale radius observed and seen in CDM simulations. For 
this reason, the inferred DM scale radius for both Abell 
383 and MS2137-23 lay at the boundary of our assumed 
prior range. Future work will benefit from weak lensing 
data, along with galaxy kinematics and X-ray data of the 
hot ICM which can each probe out to large clustercentric 
radii. 

Although not the focus of the current work, pinning 
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down the correct DM scale radius will be crucial for con- 
straining other DM mass parameters. For instance, there 
is a well-kno wn degeneracy betw e en r^e and the inner 
slope P (e.g. iGavazzi et alt (|2003f) : iGavazzil ()2005f) ). To 
briefly explore this, we have reran our analysis (for the 
coarse positioning lensing case) for both clusters with 
a Tsc of 50 and 400 kpc - factors of two beyond our 
chosen Vsc prior. We show our confidence contours in 
Figure [5l which are noteworthy. For example in the 
case of MS2137-23, if we fix rsc—50 kpc, then the best- 
fitting /3 = 0.05. However, if rsc^AOO kpc then f3 = 0.7, 
more in accordance with simulations. Interestingly, the 
rsc=400kpc scenario returns a better overall ^ 26 
than any model with rsc=100-200 kpc - even though a 
of 400 kpc is clearly ruled our by extant weak lensing 
observations. None of the other rsc=50,400 kpc scenarios 
produced values that were comparable to those seen 
with rsc=100-200 kpc. Any further knowledge of the DM 
scale radius would aid greatly in constraining f3 and de- 
termining the overall goodness of fit of the generalized 
NFW DM profile to the cluster data. 

X-ray studies assuming hydrostatic equilibrium 
(lAllen et all [200lt iSchmidt fc AIIctI [200^ and a com- 
bined strong and weak lensing analysis ()Gavazzill2005f) 
have presented data on MS2137-23 to radii much larger 
than that probed in this study. To check that the mass 
model derived from data within ~ 100 kpc do not lead 
to results at varia nce with publish ed data at larger radii, 
we have taken the IGavazzil (|2005( l results and compared 
their derived mass at large radii with an extrapolation 
of our mass models. 

Examining Figure 3 from IGavazzi (pOOl we estimate 
that from his weak lensing analysis a 2D projected mass 
enclosed between 1.6 x lO" and 1.1 x IO^^Mq at - 1.08 
Mpc using the cosmology adopted in this paper. Corre- 
spondingly, if we take all of the Ax^ < 1.0 models using 
our analysis method (the coarse positional accuracy case 
was use) and calculate the expected 2D projected mass 
enclosed at 1.08 Mpc we find values between 6.9 x 10^^ 
and 8.4 x lO^^M©, well within the expected range. 

Note that no attempt was made to extrapolate the 
mass profiles derived in our analysis to larger radii than 
the data in this paper allow, although we are acquring 
weak lensing data for a large sample of galaxy clusters 
to perform a more extensive analysis. The purpose of 
this consistency check is to only ensure that the masses 
we derive for such large radii are not too discrepant with 
existing analyses. The consistency check is satisfied and 
lends some credence to the models. 

5.2.3. Dark matter baryons interactions and Triaxiality 

The central regions of DM halos can be strongly af- 
fected by the gravitational interaction with baryons dur- 
ing halo formation. If stars form and condense much 
earlier than the DM, it is expected that the baryons will 
adiabatically compress the DM result ing in a halo that 
is steeper than that of t he original (jBlumenthal et all 
119861 iGnedin et aT]|2004f ). Ahernatively, dark matter 
heating through dynamical friction with cluster galax- 
ies can counteract adiabat ic contraction, leading to a 
shallo wer DM profile (El-Z ant et all 12004 iNipoti et all 
120031 ). The present work takes into account neither of 
the above scenarios, and if any baryon-DM interaction 
greatly changes the cluster density profile, our assumed 



param eterized gNFW profile ma y be inappropriate. Re- 
cently, IZappacosta et al.l (|2006f ) have used X-ray mass 
measurements in the cluster Abell 2589 to conclude that 
processes in galaxy cluster formation serve to counteract 
adiabatic contraction in the cluster environment. Cer- 
tainly, more observational work is needed to understand 
the interplay between baryons and DM in clusters, and 
extended velocity dispersion profiles of BCGs in conjunc- 
tion with other mass tracers at larger radii could serve 
as the best testing ground for the interplay of dark and 
luminous matter. 

Not only is there likely significant interplay between 
baryons and DM in the central regions of clusters, but 
real galaxy clusters are certainly triaxial and, if ig- 
nored, this may lead to biased parameter estimations 
and discrepancies when combining mass measurement 
techniques that are a combination of two- and three- 
dimensional. Several recent studies have considered 
the effects of halo triaxiality on observations. Using 
an N-body hydrodynamical simulation of a disk galaxy 
and performin g a 'long slit' rotation curve observation, 
iHavashi et al.l (2004) found that orientation and triax- 
ial effects can mistake a cuspy DM profile for one that 
has a constant den sity core. At the galaxy cluster scale, 
IClowe et al.l ()2004D performed mock weak lensing obser- 
vations of simulated galaxy clusters and found that the 
NFW concentration parameter recovered was correlated 
with the 3D galaxy cluster orientation. In order to in- 
vestigate the recent rash of galaxy clusters with observed 
high concentration par ameters in seeming contradiction 
to the CDM paradigm |Kneib et al. 2 003: Gavazzi 200^; 
iBroadhurst et al.ll2005bl ). IOguri et^ (|2005f ) used strong 
and weak lensing data in Abell 1689 along with a set of 
models that included halo triaxiality and projection ef- 
fects. Again, it was seen that halo shape causes a bias 
in mass (and mass profile) determination, although it 
should be kept in mind that measurements of concentra- 
tion are extremely difficult (e.g. Ha l kola e t al. 2006), and 
the recent study of iLimousin et al.l (|2006l ) has seemed to 
clear up the concentration parameter controversy for at 
least Abell 1689. 

In terms of the current work, IGavazzil ([2005) has 
pointed out that the inability of his lensing model to fit 
the MS2137-23 BCG velocity dispersion profile may be 
due to halo triaxiality or projected mass along the line 
of sight (which would increase the mass measured in the 
lensing analysis but would not b e seen in the stellar veloc- 
ity dispersion). IGavazzil (j2005f ) showed that an idealized 
prolate halo with an axis ratio of ~ 0.4 could explain the 
velocity dispersion profile in MS2137-23. Halo triaxiality 
could also explain the high concentration previously seen 
in this cluster. Again, the gap between simulations and 
observations may be bridged with respect to triaxiality 
if further steps were taken to compare the two directly. 
One step in this direction would be the publication of 
detailed density profiles for the simulations (in 3-D or 
along numerous projected sight-lines). 

The most recent DM only simulations have indi- 
cated that the standar d NFW profi l e repr esentation of 
a DM profile (and its iMoore et al.l ()1999D counterpart 
with an inner slope P ~ 1.5) can be significantly im- 
proved by slightly altering the model to a profile with 
a slope th at becomes sy s temat ically sh allower at small 
radii (e.g. iNavarro et aP ()2004l ). but see iDiemand et all 
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(|2005f )). While we have adopted the traditional gen- 
eralized NFW profile in this study, future work with 
parameterized models should move towards the latest 
fitting functions along with an implementation of adi- 
abatic contractio n as h as already been attempted by 
Zappacosta et a.1.1 (I2006D . Note, howeve r, that both 
Navarro et al.l (|2004D and iDiemand et"al1 (|2004f) have 
stated that all fitted functions have their weaknesses 
when describing complicated N-body simulations and 
that when possible simulations and observations should 
be compared directly. 

6. SUMMARY & FUTURE WORK 

We have performed a joint gravitational Icnsing and 
dynamical analysis in the inner regions of the galaxy clus- 
ters Abell 383 and MS2137-23 in order to separate lumi- 
nous baryonic from dark matter in the cluster core. To 
achieve this, we implemented a new 2D pseudo-elliptical 
generalized NFW mass model in an updated version of 
the LENSTOOL software package. This refinement is a 
natural progression from our ear Uer attempts to measur e 
the dark matter density profile (|Sand et al.ll2002ll2004D . 

For the study, we adopted an observationally moti- 
vated scale radius prior of r^c = 100 — 200 kpc. With 
strong lensing alone, we find that a range of mass pa- 
rameters and DM inner slopes are compatible with the 
multiple image data, including those with /3 > 1 as seen 
in CDM simulations. However, including the BCG kine- 
matic constraints for both systems, the acceptable pa- 
rameter ranges shrink significantly. 

We can summarize the results for the two clusters as 
follows: 

1. For the cluster Abell 383 we have found satisfac- 
tory BCG -I- generalized NFW cluster DM models 
only for our coarse lensing positional accuracy sce- 
nario. Assuming that this is reflective of the un- 
derlying cluster DM distribution, the dark matter 
inner slope is found to be /3 = 0.45°q 25, support- 
ing our earlier contention that some clusters have 
inner DM profiles flatter than those predicted in 
numerical simulations. 

2. For MS2 137-23 our model is unable to reproduce 
the observed BCG velocity dispersion profile and 
the range of accepted inner slopes therefore de- 
pends sensitively on the adopted uncertainties in 
the mass model. This may suggest an unknown 
systematic uncertainty in our analysis or that we 
have adopted an inappropriate mass model. We ex- 
plore the former in considerabl e detail, extendiii g 
the quite extensive discussion of lSand et al.l (|2004j ) . 



However, no obvious cause can be found. If, as 
we suspect, the cause lies with our adopted mass 
model, it points to the need for further work con- 
cerning the distribution of dark matter in the cen- 
tral regions of galaxy clusters. 

Future modeling efforts should include the effects of tri- 
axiality and the inffuence of baryons on dark matter. It 
is also critical to obtain high S/N extended velocity dis- 
persion measurements of more BCGs out to larger radii 
so that, in conjunction with other mass measurement 
techniques, the interplay of baryons and dark matter in 
cluster cores can be studied with a real sample. Some 
other future directions are straightforward. For exam- 
ple, the deep multib and ACS imaging now b eing done 
with galaxy clusters fBroadhurs t et al]|2005aD allow for 
literally hundreds of multiple images to be found, signifi- 
cantly increasing the number of co nstraints and allow ing 
for nonparametric mass modeling (jDiego et al.ll2004h - a 
crucial addition in case the currently used parameterized 
models do not correspond to reality. We are eager to 
find ways to more directly compare simulations with ob- 
servations so that clearer conclusions can be drawn over 
whether or not simulations and observations are compat- 
ible. This may involve measuring other properties of the 
dark matter halo rather than a sole emphasis on the inner 
slope, such as the concentration parameter, c. Simulated 
observations of nu merical simulations, such as that pre- 
sented recently bv iMeneghetti et al] (|2005| ). offer a clear 
way forward in understanding the systematics involved 
in observational techniques and the kinds of observations 
required to test the current paradigm for structure for- 
mation. 
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APPENDIX 



A GENERALIZED NFW IMPLEMENTATION IN LENSTOOL 

Here we briefly discuss the implementation of the pseudo-ellipt i cal generalized N FW proflle into the LENSTOOL 
software package. The interested reader is referred to iKneibI ()1993D : ISmith et al] (|2005i) fo r further details about the 
LENSTOOL software. Some of what follows has been presented bv iGolse fc KneibI ( 20021 ). but is reviewed here for 
continuity and clarity. 

Throughout this section we are using the thin- lens approximation with = + and x = (xi,X2) = R/rsc- 
By introducing ellipticity into the potential rather than the surface mass density we make the lensing calculations 
more tractable given that the deflection angle is just the gradient of the scaled lensing potential. Using the following 
coordinate substitution of x by Xe, 
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= y/ai xi 
— arctan (x2/a;i) 



2 , . ^2 (Al) 



where Oi and 02 are two elliptical parameters. We can calculate the elliptical deflection angle: 

de{x) = 



a{x^)y/a{coii(l)^ 

0x1 

\ ^ — = a{x^)y^a2sm(p^ 
\ 0x2 



(A2) 



The above expression holds for any definition of ai and 02 which we choose to be: 

= (A3) 

02 = 1 - e ' 

While this choice of ai and 02 do not correspond directly to the cUipticity of the potential (see lMeneghetti et al.ll200"3l 
who use a different parameterization), it does le ad to s imple expressions for standard lensing quantities, such as the 
surface mass density and shear (see Eqs. 17-19 of iGolse fc Kneih (,2002,) ) . The standard ellipticity of the potential e^p 
is related to e by 



e. = l-VT^- (A4) 

Using the standard lensing functions (see e.g. iMiralda-Escudd (|1991f) ) for the deflection angle (a), convergence (k), 
and shear (7), along with Eqn lA2l above, the projected mass density Ee(x) for our pseudo-elliptical implementation is 
simply: 

E,(f) = E{xe) + ecos2(j),{E{x,) - E(xe)). (A5) 
Likewise, the 3D pseudo-elliptical density profile can be similarly derived to be 

Pe{x, z) — p{xe, z) + ecos2(j)f{— / xp{x, z)dx — p{xe, z)) (A6) 

Jo 

where z is the direction along the line of sight. Since we are not making an effort to probe the triaxiality of our galaxy 
clusters, we will plot projected quantities whenever possible. 
Although Eqns IA6I and IA5I are general, we are working with the generalized NEW density profile 

^'^'^ " {r/rscnilX/rscW-P ^^^^ 

and the resulting surface density profile 



.7T/2 

IlgNFwiR) ^'ipcTscScX^-^ d0sme{smd + xf-^ (A8) 





when modeling the cluster dark matter halos. 

In Figure lA.ll the extent to which our implementation does not produce surface density profiles with true elliptical 
isocontours is illustrated. As the parameter e increases the surface density isocontours become more boxy and peanut 
shaped. However, at relatively low e, the isocontours are very nearly elliptical. We discuss in the following section 
to what extent our generalized NEW pseudo-elliptical mass model is an adequate description of an elliptical mass 
distribution. 

Unlike the NEW profile, the surface mass density and deflection angle of the generalized NEW profile cannot be 
calculated analytically. This greatly slows any lensing computation, especially when we need to calculate values 
over large parameter hypercubes. To limit the computing time necessary, we created a look-up table for all of the 
necessary integrals from which we interpolate when performing our lensing calculations. 

Limitations of the Pseudo- Elliptical Treament 

In this section we quantitatively investigate the range of e for which the generalized NEW pseudo-elliptical mass 
model is an adequate description of an elliptical mass distribution. As can be seen in Eigure IAtI our pseudo-elliptical 
representation can depart strongly from a true elliptical model at high ellipticities. To what degree can we consider 
out treatment of ellipticity an accurate one for representing elliptical surfac e dens ity distributions? To answer this, 
we reapply several of the quantitative measures presented by iGolse fc KneibI ()2002l ) to our generalized NEW pseudo- 
elliptical model. To get a feel for the relation between e and the ellipticity in the surface mass density, es, we plot 
several values in Eig. [T] 
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We quantify the degree of boxiness by measuring the distance, d r, between a su r face d ensity contour and a real 
elHpse with the same major and minor axis radii (as was done in iCxolse fc Knelbl (I200I ): see their Figure 6 for a 
geometrical illustration of Sr). In Fig. [5] we plot 6r/r as a function of e for several values of the inner slope, /3, and a 
variety of r/rsc- If we desire our pseudo-elliptical implementation to be within 10% of a true elliptical surface density 
distribution for r/rsc < 10, then values of e < 0.2 arc appropriate, especially for DM halos with steep inner slopes. 

One unphysical consequence of introducing ellipticity into the potential is that the surface mass density can become 
negative, especially near the minor axis where cos{2(j)^) = —1. In Figure P^. II we plot the distance along the minor axis 
at which becomes negative for several inner slopes. If we wish to restrict ourselves to values of e where the surface 
density does not go negative for r/r^c < 10, then we must restrict ourselves to values of e less than approximately 
0.25. 

In summary, in order to be within 10% of a true elliptical surface mass distribution and to have positive values of 
the surface mass density for r/rgc < 10 we must restrict our use of the pseudo-elliptical gNFW parameterization to 
values of e < 0.2, well within the model values of e for the clusters studied in this paper. 
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Fig. 6. — Illustration of projected density isocontours for the pseudo-elliptical generalized NFW parameterization with r/rsc = 10.0 and 
/3 = 1.0. Note that as e gets larger, the projected density contours become more dumb-bell shaped. 
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Fig. 7. — The ellipticity (talcing the minor and major axis positions and assuming ex; = 1 — b/a) of the projected density, S, as a function 
of the ellipticity in the potential for different values of /3. 




Fig. 8. — Sr/r as a function of e for a variety of pseudo-eUiptical generalized NFW models with different inner slopes, /3. This simply 
characterizes the deviation of the projected density from an elliptical model for various r/rsc- 
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Fig. 9. — Distance from ellipse center along the minor axis at which becomes negative. Several example values for different inner dark 
matter density slopes, (3, are plotted. 



